FAX NO. b208283384 14 


maw nr? nn r-n i 

rrci 


< n . 

ID- 


OlJ'T'A A T f 

Sniu bi 


6 S"0- Cc '/- ^ 


X f-A f p !' C tj I— 

Cl 

Final Technical Report: Airborne Astronomy Program 
Grant NAG 2-1074 

The Circumxtdlar Environment of Low Mas# Star Forming Regions 
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1- Introduction 


Since the 1980*$, our understanding about the intcr-rclationship between the col- 
lapsing cloud envelope and the disk has been greatly altered While the dominant 
star formation models invoke free fall collapse and a density profile, other star 

formation models are possible. These models invoke either different cloud starting 
conditions or the mediating effects of magnetic fields to alter the cloud geometry 
during Collapse. To test these models, it is necessary to understand the envelope‘s 
physical structure. The discovery of disks, based on millimeter observations around 
young stellar objects, however makes a simple interpretation of the emission com- 
plicated. Depending on the wavelength, the disk or the envelope could dominate 
emission from a star. In addition, the discovery of planets around other stars ha* 
made understanding the disks in their own right quite important. Many star forma- 
tion models predict disks should form naturally as the star is forming, in many cases, 
the information we derive about disk properties depends implicitly on the assumed 
envelope properties. How to understand the two components and their interaction 
with each other is a key problem of current star formation. 


To .separate out the two components, it is necessary to observe sources at a variety 
of wavelengths, over a high spatial resolution. As mentioned above, depending on 
the physical conditions in the disk and the envelope, each will contribute differ rent 
amounts of flux at each measured wavelength. In addition, high spatial resolution 
should allow us to resolve the envelope emission. Envelope emission typically extends 
over many hundreds of AU, whereas the disk will be confined to a much smaller area. 
By combining high resolution multiwavelength data with radiative transfer models, 
one can assess the likely contribution of both the disk and the envelope. The envelope 
contribution dominates the far-infrared emission, with as much as 80% of the total flux 
corning from the cool dust found in tire outer portions of the envelope. This emission 
can be detected and mapped with high resolution observations, such as the Yerkcs 
Far-IniraRed (Flit) Camera on board the Kuiper Airborne Observatory! KAO). In 
1993 and 1995, we obtained far-infrared (60. 100, 1(50, 200 micron) images for a 
number of young stellar objects (VSOs) with the KAO prior to the KAO's shutdown. 


In the proposal covered by the grant (NAG-2-1074), we proposed to analyze the 
dat a, acquire additional data at other wavelengths for models, and to develope source 
models for the YSOs under study. Based on the results of these models, we hope to 
define the range of physical conditions; that are expected in young stellar systems of ' h/^J/‘ 
different ages and masses. Through the A we were a ?>?e 
to accomplish manv of our project goals, jm if 
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II. Observational Data Acquired for Source Models 

Far-Infrared Images: Iti KY 93 and FY95, we were awarded two flights each 

year to study the envelope arounds low mass protostars in Taurus. Working with 
the Ycrkcs FIR. camera, wc obtained data on a total of 13 sources. The Ycrkes FIR. 
camera, contained an tit.) element, array. Our standard observing method was to dither 
the array around the desired source position, and to take a series of images for later 
coadding. Initial data reduction indicated that most of the sources were compact, a„s 
would be expected for low luminosity protostars surrounded by an infalling envelope. 
Such objects will not tend to heat up a lot of the envelope, and what dust is heated 
will be concentrated towards the center. However, two sources ( hi 52 7 and L155INE) 
appeared to show extended emission. This was particularity interesting since these 
objects may be the the youngest objects in the sample. The claim of youth is based 
on their spectral energy distributions (SF.Dx), which are. peaked very strongly in the 
FIR. Such a peak implies that, the object is still deeply enshrouded in dust. 

Our initial data, reduction was based on smoothing the data with a relatively 
broad smoothing function (i.c. a Gaussian beam). The smoothing resulted in the 
data being far from the diffraction limit. To address this situation, we undertook 
( with the Yerkcs group) to re-reduce the data using improved convolution methods. 
Due to the desire, we had to spend considerable time to achieve the diffraction limit. 
These new results are still being analyzed, however, they confirm the detection of 
extended emission around LlfxilNE and L1327. The final models of the KAO results 
are still in preparation. 


To supplement the KAO data, we obtained observations at related wavelengths 
for a number of the sources in our KAO dataset. 

Mid-Infrared spectra: A large, sample of T Tauri stars was studied by Cohen 

and Witteborn in 1985 at 10 pm. To extend this survey, we obtained R=200 spectra, 
for several T Tauri stars, and a more extended sample of Herbig Ae/Be stars. These 
stars have a variety of 9.8 /mi silicate emission. This emission feature reveals that 
optically thin silicate dust is in the system. Such dust is similar (spectroscopically) 
to dust seen around some Vega-like objects like Beta Pic. 

The passage of cornet Hale Bopp provided a.n opportunity to measure the spec- 
troscopic properties of comets for comparison with our T Tauri and Herbig Ae/Be 
sample. The cornet spectra was found to have some features that were time variable 
The comet spectra reveal the presence of grain oomponenents that might be pre-solar 
(based on comparison with IDP spectra). We hope to tie the Comet Hale Bopp work 
to our FIR/M l R studies in our final models. 

Mid-Infrared images: Using the MIRAC2 (Mid-Infrared Array Camera), we 

imaged several of the embedded YSOs at 10 /mi. These observations provide high 
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sensitivity high spatial observations of the circumstellar environment on size scales 
of 1 r>0 AU. Most of the observations arc consistent with the mid-infrared emission as 
corning from a disk. Extended emission is possible, since we found several examples 
in related studies of Her big Ae/Bo stars. However, in general, the T Tauri/ embedded 
low mas* stars are too cool to heat an extended envelope unless some unusual grain 
heating mechanism is invoked . 

VVe obtained 10 and 20 /r m images of the Trapezium duster, with its proplyds, 
which we are currently analyzing to provide a. comparison with our isolated stars. 
Initial results are being published in Deutsch ot al (1999). 

SubmiHimeter/millimeter photometry: We focused on Taurus embedded 

VS OS and the Vega-like stars. In the case of the Taurus sources, which includes our 
KAO sample, we have obtained a complete survey of the 450-2000 /.mi emission using 
the JCMT lo meter telescope. For the Vega-like stars, we used the C$0 10m and 
I RAM 30m telescopes to get 300 and 1300 /j m photometry. The data will be written 
up in two papers. 

Submillimeter/Millimeter Interferometry: 

A number of our sample stars have data already available in the literature. For 
t hose sources not already observed, we obtained 1.3 mm interferometry data for the 
brightest 16 YSOs in our Taurus sample. These observations extended a previous 
survey we had done. We detected 14/16 sources (Wannicr ct al. 1999 in prep). Of 
particular interest, we discovered evidence for a couple of binaries in our sample 
(L1551NE). 'This dataset is to bo combined with published interferometer 2.7 mm 
data, reported in several surveys. Wc also were awarded time to study a. sample ot 
Taurus YSOs with the JCMT-CSO Si ngle- Beam- In ter feroi noter , and got data on DC 
Tau , among other sources. 

Summary: Wc have obtained the complete $ED from 10 microns out to 

1.3mm for all of our sources. We have the FIR imaging data, processed to reveal the 
maximum angular resolution possible, which allows us to model the disk lo model 
the disk, we have the high resolution millimeter interferometry data. 
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111. Modeling Techniques Developed for Source Models 


At the start of our KAO research, we had one primary modeling technique: a .1-0 
radiative transfer code that allowed us to explore a variety of density gradients and 
grain properties. We also developed a way to including disk emission in our model. 
With these models, we ran explore many types of disk/cnvelope models. 

We identified two other radiative transfer codes that would be useful for our 
models: 

]..) Inclusion of Polycyclic Aromatic Hydrocarbons (PAHst and small grain heat- 
ing for the models. This code, developed by Siebenmorgcn, Krugel. and Mathis 
(1992), allows us to do spatial models of the PAH emission and small grain emission. 
Such effects are critical to understanding the emission from a number of the Herbig 
Ae/Be stars we have under study. 

2.) A fast 2-D code: Developed at Jena by Dr. Men'shchikov, it can handle 
both deeply embedded sources as well as visible disks. It. has been applied by Dr. 
Men’shchikov to Li 051 and 1IL Tan. We intend to apply it to our more deeply 
embedded sources (L1551NE and L1527). 

Overall, there are three important parameters we examine: 

1) The envelope geometry and temperature structure. 


This is best done with the SED and spatially resolved data at far-irifrared/sub- 
1 1 1 i 1 1 i meter wa w lengths. 

2) The disk geometry and temperature structure. 

This is most sensitive to the mid- infrared emission, due to the role of small grains 
or PAHs. High resolution data at 10 microns provides the best constraints. 

3) The dust, properties in both the envelope and the disk. 

The models allow a- tradeoff (to a limited extent) between dust, properties and 
the inner stucrure of the disk/envelopc system. However, the presence or absence of 
silicates is strongly limited by the 10 micron spectra. 
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IV. Results to date 


Combined with our iuif Ul KAO program, we find* 

1. ) The vast majority of embedded YSOs in Taurus are compact at 100 microns. 
The models most consistent with our data and other observations are either dominated 
by disk emission (for sources like GO Tau) or envelopes that have relatively steep 
density gradients ( r — 1 -5to-2.0) < 

2. ) The subinillimetcr/millimeter photometer suggests that envelope models 
(such as those of Hartmann ft al. (1993) are very successful. However, in a few 
oases) they are clearly missing a component in their models. Disk emission plays an 
important role and must be considered when predicting the overall emission. (Butner 
cl al. 1999b(in preparation)). 

:j) In the two cases where we seem to have extended emission, wo have to investi- 
gate other possible source models than a Shu collapse. This work is still in progress. 
We will considering models such as those of Hartmann cl al. (1995), which invoke a 
flat sheet rather than spherical collapse. 


V. Future Work 


We still doing the final modeling for the revised FIR maps, and anticipate having 
two or three papers done by the end of the year. These papers will present our dal.* 
and the complete set of parameter space that we explore. In particular, we wish 
to assess the role of flat sheet collapse models (see point 3 above) in explaining our 
Taurus FIR data. 
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VI. Papers in Preparation or Published 


The following papers were supported in part by the Airborne Astronomy Grant 
NAG 2-1074. The support provided salary, travel funds, or equipment used for data 
reduction and analysis. 
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H. .J. Walker, V. T sikoudj, C. A. Clayton . T. Gcballc, D. H. Wooden, H. M. Butuer, 
“The Nature of the unusual source IRAS 18530+0817,” 1997. A&A. 323, 412. 

H. M. Butncr, H. J. Walker. D. II. Wooden, & F. C. WiLteborii. “1997." Examples of 
Comet-Like Spectra Among 0 Pic- Like Stars, in Astronomical and Biochemical 
Origins and the Search for Life in the Universe , ed. C. B. Cosmovici, S. Bowyer, 
Sc D. Wcrthimcr (Bologna: Editrice Compositori), p.,149. 

C. H. Mori arty-Schie veil Sc H. M. Butner, “A Sub-Millimeter- Wave ‘Flare’ from GG 

Tan?,” 1997. Ap. 474. 70S. 

Selected Unrefereed Articles 

A. Natta, L II. Butner, “Resolving Disks in YSOs.” 1996, in Infrared Space fntcr- 
ferometry. ed. A. Albordi, T. dc Graauw, C. Eiroa, C. J. Schalinski, and H. 
Timm son (Space Science Reviews), in press. 

Papers in Preparation 

H M. Butner. G. II. Moriarty-Schicvcn, P. Wannier, “Spectral Energy Distributions 
in Taurus, 1999, Ap.J.. in preparation. 

H. M. Butner, G. M. Moriarly- Scliieven, el. ul. . “FIR Imaging of Low' Mass Sources 
in Taurus,’ 1999, Ap.J., in preparation . 

L. K. Deiitsch, M. Kassis, N. Smith, I. Hora, G. G. Faaio, H. M. Butner, W. F. 
Hoffman. A. Duval. “Wide-Field Thermal Imaging of the Qrion Nebula at High 
Spatial Resolution." 1999, Ap.J., in preparation. 

P. G. Wannier, J. Keene. G. Moriarty-Schieven. and II. M. Butner. “13 MM inter- 
ferometry of Taurus Sources,” .1999, Ap. in preparation. 

D. II. Wooden, H. M. Butner, Sc K. C. Witteborn, “A Mid-Infrared Survey of T Tauri 

Stars.” 1999, Ap. J ., in preparation. 


FAX NO. 5208283384 14 


maw m n n rnt i n , 

iinr-u/-^ rru 100^ 


'aruu bl \t 


< 

References 

Below we include the references for articles mentioned in our report. 

Butncr, H. M.. Walker, H. J, Wooden, D. H., k Witteboro, F. C. 1907, in Astronom- 
ical and Biochemical Origin* o.nd the. Search for Life in the Universe , ed. C. B 
Cosmovici. S. Bowver, k D. Wcrthirnei (Bologna: Fditrice Compositori), p.149. 

Butner, H. M., Walker, H. J., Wooden, D. II., k Witteboro., F. C. 1999a, Ap. J ,, in 
preparation. 

Butncr, H. M., Mori art y-Schieven, G. H., and Wannier, P. 1990b. Ap.J., in prepara- 
tion. 

Cohen, M., Witteborn, F. C 1985, A pJ , 294, 345. 

Deutach, L. K.. Hora, J. L., Bu trier, H. M. , Hoffmann, W. K. k Fazio, G. G. .1.995. 
in Circtimsidlar Matter 1994 > ed- G. D. Wait, (Kluvver: Dordrecht), p.S9. 

Hartmann, L, Kenyon. S. and Calvet. N. 1993, ApJ, 407. 219. 

Hartmann, L., Cal vet, N., and Boss, A. 1996. A p.J. % 464, 387. 

Marsh, K. A- ; and Mahoney, M. J. 1992, ApJ , 395, L.115. 

Moriarty-Schieven, G. H. and Butncr, H. M. 1997, Ap. J., 474, 70S. 

Wooden, D. H., Butner, H. M., k Wittcborn, F. C. 1999., Ap. in preparation. 

Wannier, P. G.. K eerie, J., M oriarty-Schicvcn , Butner/H. M. 1999, Ap. in prepa- 
ration. 


